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5^ . In view of the forthcoming experiments at high energy e^e" col- 

1-^ \ liders, LEP2 and beyond, the process e^e~ — > W~^W~ — > 4/ has 

been investigated for off-shell VF's. After deriving the matrix element 
^H. for the process within the helicity amphtude formahsm, QED radia- 

te I tive corrections to the initial state leptons have been included at the 

leading logarithmic level in the framework of the electron structure 

functions. Cross sections and distributions are computed by means 

of both a semi-analytical and a Monte Carlo code, allowing for sev- 

C^ '_ eral cuts of experimental interest. A Monte Carlo event generator has 

also been built, for simulation purposes. Some illustrative numerical 
results are included and shortly described. 
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In the context of the experiments planned at e'^e~ coUiders such as LEP2 
and beyond, the process e+e~ — >• W^W~ ^ 4: f plays a special role since 
it will be used for an accurate measurement of the ly-mass and hopefully to 
test the non-Abelian trilinear couplings of the standard model as well as the 
effects of non-standard "anomalous" couplings between the gauge bosons. 
In order to achieve these experimental goals, accurate predictions for the 
observables of the reaction e^e~ —>■ W~^W~ -^ 4: f are mandatory and in 
particular the contribution of all radiative corrections should be kept well 
under control. 

For the above reasons a number of studies devoted to four fermions pro- 
duction in electron-positron collisions, including the main effects of radiative 
corrections, have recently appeared in the literature and they can be roughly 
grouped in two main lines, which are in some sense complementary: 

• semi-analytical approaches 0, |^, 

• Monte Carlo approaches ||^, ^, |^ . 

Previous calculations of the higher order QED corrections to e+e^ -^ W^W~ , 
including the influence of the finite W width, were performed in [^ ^ . The 
complete study of the reaction e~^e~ -^ 4 fermions clearly represents a multi- 
step program and the aim of this letter is to summarize our novel theoret- 
ical approach to it, which includes the description of e^e~ -^ W^W~ — > 
4 fermions with initial state electromagnetic radiation in the leading log ap- 
proximation. A complete analysis of the background processes and of the 
exact 0{a) radiative corrections are currently under investigation. Our main 
goal in this letter will be to present an illustrative sample of the numerical 
results obtained so far by means of a semi-analytical and a Monte Carlo 
program. We start with a brief description of the basic ingredients. 

We have undertaken a general study of the process e'^e^ -^ 4 fermions 
with the aim of producing and discussing a wide number of realistic distri- 
butions. As a first step we briefly introduce some basic concepts and apply 
them to the double-resonating diagrams e^e~ -^ W^W~ -^ 4 fermions. The 
main emphasis in our approach is on a covariant description of the kinemat- 
ics, which allows us to translate any sort of kinematical cut directly into an 
analytical constraint on the integration variables. As for any 2^4 process 
we can introduce 15 invariants Sij = —{eipi + (-jPjY where pi^i = 1, ... ,6 



are all the momenta assumed to flow inwards and Cj = ±1 for incoming 
(outgoing) particles. There are 6 linear and 1 non-linear constraints among 
the Sij with a total of 8 independent invariants, one of which is of course the 
invariant mass of the event, s. Using standard techniques all the components 
of the momenta (up to a sign ambiguity) can be expressed in terms of the Sij 
and thus in terms of the independent invariants. This is why we can express 
any kinematical cut as an analytical constraint on the integration variables 
and why our mult i- dimensional phase space can be treated exactly. 

For e+(j9i)e^(p2) -^ fiiPs) + ■ ■ ■ fiiPe) the independent variables used 
to describe the process are, besides s = — (pi + ^2)^, the following ones: 
zi = xi3 + X23, Z3 = xi5 + 0:25, /i? = 3:34, fil = X5e, l4 = ^35, n = 0:13 and 
T2 = a;i4 where we have defined XijS = Sij. In the cm. system they have a 
rather simple interpretation, namely 

^ = 1,3, (1) 

where Ei are the energies of particles 1 and 3, which correspond to the down 
fermion coming from the decay of W~ and to the up fermion coming from 
the decay of W~^, respectively; 

f^-=-, ^ = 1,2, (2) 

s 

where Si and S2 are the invariant masses of W^ and W~ respectively; 

r, = ^(l-cos^+), ^ = 1,2, (3) 

where i^f is the scattering angle of the i-th particle with respect to the 
incoming positron beam; 

/^3 = —, (4) 

where S13 is the invariant mass between the particles 1 and 3. Given the in- 
dependent variables, the further choice of the solution of the non-linear con- 
straint allows the reconstruction of the kinematics in the laboratory frame. 
The explicit solution of the kinematics allows us to translate several cuts 
of experimental interest in terms of boundaries of the independent variables 



phase space; such a translation is of great importance from the point of view 
of numerical integration when adopting a semi-analytical approach. 

When including initial state radiation, two more independent variables 
Xi^2; which represent the energy fraction carried by the incoming electron 
and positron after the radiation, have to be added to the previous set. In 
this case the above independent variables are referred to the cm. system 
of the hard scattering process and in this frame the dependent variables are 
derived by explicit solution of the kinematics. By means of a Lorentz boost 
the whole kinematics in the laboratory system can finally be reconstructed 
thus allowing the boundaries of the phase space in the presence of realistic 
cuts to be defined also for the radiative process. 

Our matrix element contains for the moment only the three double- 
resonant diagrams in Born approximation, but there is no limitation of the 
method that will prevent us from a straightforward generalization to all the 
background processes. The strategy adopted is to use one of the many for- 
mulations of the so-called helicity method 0], which has the advantage of 
working directly with invariants. Roughly speaking, each diagram can be 
written as 

di{X,pf,p'f) = {2nfig^r]iKiAiV^x^inUxUp^TfV-p^Up^,Tf>V-p^, (5) 

where A, pf and pf are helicity labels, rji is the sign of the diagram, Kj stands 
for a collection of constants and Aj for the collection of propagators. Finally 
Tin, Tf and Tf are strings of 7 matrices. Once each diagram is rewritten 
as a trace, we explicitly insert for the bilinears the corresponding expression 
(valid up to an overall phase), as for instance 

Mp)M<i) = -Ij^,{^-^^')^^ (6) 

and the trace operation will convert everything in terms of invariants, which 
in turn are related to the integration variables. According to the above 
procedure the three diagrams of the process e~^e~ -^ W'^W~ ^ 4: f have been 
computed by means of SCHOONSCHIP H and FORTRAN-coded to get the 
squared matrix element numerically . The effect of the electromagnetic initial 
state radiation is treated in the language of the QED structure functions [jlO|, 



11 , where the corrected cross section can be written in the following form: 



a{s) = / dxi dx2 D{xi, s)D{x2, s)d[PS] . (7) 

In the latter equation (i[PS'] denotes the volume element in the 7- dimensional 
phase space connected to the independent variables of the hard scattering 
reaction in the cm. frame; D{x,Q'^) provides the probability of finding, 
inside a parent electron, an electron of energy fraction x with virtualness 
Q^. The explicit expression which we adopt to treat initial state radiation is 
given by 



32 ' 



r(i + i/?) 
1 



^(^'^) = Z^TT^^VT^^ ^(^ - ^)"'' - ji^ + ^) (8) 



4(1 + x) ln(l - x) + 3(1 + x) Inx - 4 5- x 

1 — X 



with 



P = 2-{L-1), V = 2-L, (9) 

TT 71 

where L = In {s/m?) is the coUinear logarithm, 7^; is the Euler constant and 
r(z) is the gamma function. The first exponentiated term of Gribov-Lipatov 
form describes multiphoton soft radiation, the second and third ones hard 
collinear bremsstrahlung. In the structure function the arbitrary scale Q^ is 
set equal to s and the hard photon terms are taken with the coefficient 77 
instead of (3 in order to treat the hard radiation in leading log approxima- 
tion (in this approximation, problems connected with gauge invariance are 
avoided) . 

As a consequence of the high dimensionality of the integration, a Monte 
Carlo strategy could be in principle competitive with a deterministic numeri- 
cal procedure, as far as the balance between number of calls and accuracy are 
concerned. Both of the approaches have been adopted to cross-check and put 
under control the numerical results. The Monte Carlo integration (based on 



the random numbers generator RANLUX |T^) has been performed according to 



the importance sampling technique in order to cure the peaking behaviour of 
the integrand. This strategy has been helpful also for the semi-analytical ap- 
proach, which has been carried out by means of the NAG [O routine DOIGDF. 



As a result, the deterministic procedure turns out to be more efficient with 
regard to a totally inclusive set-up and in the presence of cuts on the final 
state invariant masses and/or energies and/or scattering angles of the out- 
going fermions. On the other hand, the Monte Carlo allows for a complete 
and more fiexible control of any kind of experimental cut. Furthermore, the 
Monte Carlo routine can also be used as an event generator for simulation 
purposes. In this case the events, defined as the components of the four final 
state particles momenta, plus the radiative variables xi and X2, plus ^/s, are 
stored into proper n-tuples. A detailed description of the formulation and of 
the technical solutions adopted will be presented elsewhere. 

In the following we show some illustrative results obtained by means of 
the two numerical procedures described above. Our set of input parame- 
ters is Mz = 91.1887 GeV, Mw = 80.22 GeV, Tz = 2.4974 GeV, Tw = 



2.08 GeV PI 



In Fig. 1 the total cross section, for W decaying in all possible final states 
(excluding top production), in Born approximation (dashed line) is com- 
pared with the initial state QED-corrected one (solid line) as a function of 
the cm. energy for a fully extrapolated set-up. Concerning the tree level 
prediction, the dashed line result has been obtained by a two-fold integration 



of the formula given in ref. ||T^ while the solid line derives from a nine-fold 
integration of the e~^e~ -^ W^W~ -^ 4/ matrix element performed by our 
semi-analytical code. The Monte Carlo results are represented by the open 
circles, showing a fully satisfactory agreement in both the Born and QED- 
corrected case. For such a configuration, the effect of initial state radiation is 
to reduce the peak cross section by about 10%, to slightly shift the peak po- 
sition towards higher energies and to determine the appearance of a radiative 
tail as a typical convolution effect. For the QED-corrected total cross section 
we compared our results, obtained by properly adapted structure functions, 
with the values quoted in [Q and found very good agreement. 

Figure 2 shows the W~ invariant mass distribution da/dM^ (M^ = 
— (gi + ^2)^) at the Born level for two cm. energies in the LEP2 regime, 
namely y/s = 175 GeV (dotted line) and ^/s = 190 GeV (dashed line). In or- 
der to single out the contribution of the initial state radiation, the corrected 
invariant mass distribution at ^/s = 175 GeV (solid line) is also plotted for 
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comparison. As in Fig. 1, the Monte Carlo predictions given by the markers 
fully agree with the semi-analytical results used to produce the continuous 
lines. 

The effect of some realistic experimental cuts is shown for the total cross 
section (Fig. 3) and the invariant mass distribution at ^/s = 175 GeV 
(Fig. 4) of e^e~ — > W~^W~ — > 4/ where both W^s are required to decay 
hadronically. All the curves include the effect of the initial state leading 
log QED corrections. The experimental conditions considered are the fol- 
lowing: no cut (solid line), cuts on the energies of all the four final state 
fermions -E^m — 25 GeV (dashed line), constraint on the invariant masses 
of both W^s 77 < M± < 83 GeV for the total cross section and on the W~^ 
invariant mass 77 < M+ < 83 GeV for the W^ invariant mass distribution 
(close-dotted line), acceptance cuts on the scattering angles of the produced 
fermions 40° < 'di < 140° (wide-dotted line). Also for such realistic configu- 
rations the agreement between the semi-analytical (lines) and Monte Carlo 
(markers) integrations is fully satisfactory pointing out the flexibility of our 
semi-analytical code. 

In Figs. 5-9 we show some distributions of experimental interest for LEP2 
physics, which have been obtained processing the n-tuple created by the event 
generator at ^/s = 175 GeV taking into account the emission of initial state 
radiation. They refer to a sample of 10^ four fermion events and show the 
effects on the considered distributions of three different selection criteria: 
fully extrapolated set-up (white histogram), cuts on the scattering angles of 
1, 2 particles coming from the decay of the W~ boson (40° < -(91,2 < 140°) 
(dark histogram) and cut on the invariant mass of the W~ boson (77 < 
M_ < 83 GeV) (grey histogram). The effect of the cuts is shown for the W~ 
angular distribution (Fig. 5), the energy distribution of the fermion from W~ 
decay (Fig. 6), the distribution of the relative angle of the decay products of 
the W^ boson (Fig. 7), the scattering angle of the fermion produced by W~ 
decay (Fig. 8) and, finally, the photon energy distribution (Fig. 9). 
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Figure Captions 



Figure 1. The total cross section of e^e^ -^ W^W~ -^ 4/ without (dashed 
hne) and with (sohd hne) initial state leading log QED corrections as a 
function of the cm. energy and for W^s decaying into all possible final states 
(top excluded). The open circles are the results of our Monte Carlo program 
in comparison with the integration of the analytic formula of ref. [15] for the 
tree level approximation and with the predictions of our semi-analytical code 
for the radiative case. 

Figure 2. The W~ invariant mass distribution da/dM^ (M^ = — (gi + ^'2)^) 
at the Born level for two cm. energies ^/s =175 GeV (dotted line) and 
^/s = 190 GeV (dashed line). The solid line is the QED-corrected invariant 
mass distribution at ^/s = 175 GeV. The Monte Carlo integration results are 
given by the markers. 

Figure 3. The effect of some experimental cuts on the total cross section of 
e"'"e~ — >■ W^W~ -^ 4/ including initial state leading log QED corrections as 
a function of the cm. energy, and for both W bosons decaying hadronically. 
The Monte Carlo integration results are given by the markers. 

Figure 4. The same as in Fig. 3 for the invariant mass distribution at ^/s = 
175 GeV. 

Figure 5. The effect of specific selection criteria on the W~ angular distri- 
bution at ^/s = 175 GeV with QED corrections included: no cuts (white 
histogram), cuts on the scattering angles of 1,2 particles coming from the 
decay of the W~ boson (40° < 'i9i,2 < 140°) (dark histogram) and cut on the 
invariant mass of the W~ boson (77 < M_ < 83 GeV) (grey histogram). 

Figure 6. The same as Fig. 5 for the fermion energy distribution. 

Figure 7. The same as Fig. 5 for the relative angle of the decay products of 
the W~ boson. 

Figure 8. The same as Fig. 5 for the fermion scattering angle. 

Figure 9. The same as Fig. 5 for the photon energy distribution [E^ < 
25 GeV). 
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